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We present an analysis of the MSW neutrino oscillation solutions of the solar neutrino problem in the 
framework of two-neutrino mixing in the enlarged parameter space (Am 2 , tan 2 9) with 9 £ (0, ■§). 
Recently, it was pointed out that the allowed region of parameters from a fit to the measured total 
rates can extend to values 9 > j (the so called "dark side") when higher confidence levels are 
allowed. The purpose of this letter is to reanalize the problem including all the solar neutrino data 
available, to discuss the dependence on the statistical criteria in the determination of the CL of the 
"dark side" and to extract the corresponding limits on the largest mixing allowed by the data. Our 
results show that when the Super-Kamiokande data on the zenith angle distribution of events and 
the spectrum information is included, the regions extend more into the second octant. 
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In a recent work, de Gouvea et al. [1] have pointed out 
that the study of two-active neutrino oscillations in the 
framework of the MSW [2] solutions to the solar neu- 
trino deficit as done traditionally on the (Am 2 , sin 2 (29) 
parameter space is incomplete since it covers only the 
range < 9 < ~ (which they denote as "light side"). 
By fitting the data on the total rates meassured at Chlo- 
rine [3], Gallium [4,5] and Super-Kamiokande [6] experi- 
ments, they claim that the allowed region of parameters 
extends to the region with 9 > ? (denoted as the "dark 
side") at some reasonable confidence level (CL). In fact, 
to our knowledge, the need of extending the mixing pa- 
rameter space was first discussed in Ref. [7] where the 
mixing variable tan 2 9 was introduced to chart the full 
mixing range < 9 < ^. 

In this note we revisit the problem of the extension 
of the "dark side" , that we will denote simply as second 
octant, of the solar neutrino parameter space after in- 
cluding the effect of all other solar neutrino observables. 
In particular in this analysis we use all measured total 
event rates as well as the 825-day Super-Kamiokande 
data on the zenith angle dependence and the recoil elec- 
tron energy spectrum of the events. We also discuss the 
dependence on the statistical criteria used in the con- 
struction of the allowed regions. Our results are sum- 
marized in Tables I and II where we show the maximum 
allowed values of tan 2 9 at the 90 and 99% CL when the 
different observables are included as well as the CL for 
the second octant of the parameter space. We show that 
when adding the zenith angle and the spectrum infor- 
mation in the analysis the regions extend more into the 
second octant. In particular for the LMA solution, this 
behaviour is mainly driven by the data on the zenith an- 
gular dependence since the best fit point for the zenith 



angle distribution is in the second octant as pointed out 
in Ref. [8] . For details on the data and statistical analysis 
employed in this study we refer to our detailed work on 
Ref. [9] and references therein. 

Let's first recall that in the framework of two massive 
neutrinos, the weak eigenstates (v e and v x for the solar 
neutrino problem) can be writen as a linear combination 
of the mass eigenstates v\ and vi 



v e = cos v v\ + sin tt v-2 
v x = — sin 6 V\ + cos 9 V2 , 



(1) 



where 9 is refered to as the mixing angle in vacuum. 
The mass-squared difference is defined as Am 2 = m| — 
to 2 . For the solar neutrino problem v x can label either 
an active neutrino x = fi, r or an sterile neutrino. In 
what follows we restrict our discussion to oscillations into 
active neutrinos. For oscillations into sterile neutrinos 
large mixing solutions are not allowed [9]. 

By inspection of the symmetry properties of Eq. (1) 
one sees that the full parameter space can be exhausted 
by using the mass-squared difference Am 2 always pos- 
itive and the mixing angle in the interval < 9 < 5. 
In the case of vacuum oscillations, moreover, the transi- 
tion probabilities as directly derived from Eq. (1) can be 
writen in terms of sin 2 (29) and therefore they are sym- 
metric under the change 9 — > 9 — ^ so each allowed value 
of sin 2 (29) corresponds to two allowed values of 9. 

On the other hand, in the case of the MSW solutions 
the transition probability takes the form 



r> T^Sun jjEarth 

-fee — -r e l r \e 



pSun r>Earth 
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(2) 



where Pf™ is the probability that a solar neutrino, 
that is created as v e , leaves the Sun as a mass eigen- 
state fi, and p^ e arth is the probability that a neutrino 
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which enters the Earth as v\ arrives at the detector as 
u e [10]. Similar definitions apply to P e s 2 u ™ and Pg arth . 
For P^ arth we integrate numerically the evolution equa- 
tion in matter using the Earth density profile given in 
the Preliminary Reference Earth Model (PREM) [11]. 

The quantity P^ un is given, after discarding the oscil- 
lation terms, as 



Sun 
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Sun 



\ + {\-PLz)cos(20 mfl ) (3) 



where Plz denotes the standard Landau- Zener probabil- 
ity [12] and # mj o is the mixing angle in matter at the 
neutrino production point: 
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with A = 2V2G F EN e (r ) where N e (r Q ) is the num- 
ber density of electrons in the production point, E is the 
neutrino energy, Gf is the Fermi constant and r res is the 
resonant point (Am 2 cos 29 — A res ). This probability is 
clearly not invariant under the change 9 — > -| — 9 as res- 
onant transitions are only possible for values of 9 smaller 
than j as seen in Eq. (4) what considerably supresses the 
transitions for 6 > j. For this reason most of the ear- 
lier papers on the MSW effect that considered the two- 
neutrino mixing case used (Am 2 , sin 2 29) as parameters 
in the fitting procedure [13] with 9 € [0, j]. In some cases 

the parameter space was represented as , <fn [14,15] 

COS (/ 

but still assuming 9 to be in the first octant. 

Due to the fact that this choice does not exhaust the 
full parameter space once matter effects are included, 
other representations have been used to show the en- 
larged space particularly in the framework of three- 
neutrino [7,16] and four-neutrino oscillations [8]. In this 
way two suggestions to parametrize the mixing angle 
have been made in the literature: sin 2 9 and tan 2 9 with 
9 € [0, §]. For the sake of comparison with the analy- 
sis in Ref. [1] we choose to show our results in terms of 
tan 2 9. 

We now turn to our results. We first determine the 
allowed range of oscillation parameters using only the 
total event rates of the Chlorine, Gallium and Super- 
Kamiokande experiments. We have not included in our 
analysis the Kamiokande data [17] as it is well in agree- 
ment with the results from the Super-Kamiokande exper- 
iment and the precision of this last one is much higher. 
For the Gallium experiments we have used the weighted 
average of the results from GALLEX [4] and SAGE [5] 
detectors. Our results are shown in Fig. 1 where we 
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FIG. 1. Allowed regions in Am and tan 8 from the mea- 
surements of the total event rates at Chlorine, Gallium and 
Super-Kamiokande (825-day data sample) in the two differ- 
ent statistical criteria. The darker (lighter) areas indicate the 
99% (90%)CL regions. The best-fit point used to defined the 
regions are indicated by a star. 



show the allowed regions. We choose to show our al- 
lowed regions for the "conventional" 90 and 99 % CL. 
It is obvious that increasing the allowed CL would lead 
to larger regions. The CL for the second octant and the 
maximum values of tan 2 9 for which the LMA and LOW 
solutions are allowed at those CL can be found in Ta- 
bles I and II. Furthermore, in Fig. 2 we show the value 
of X 2 rli „(tan 2 6<) - x 2 min m the LMA and LOW regions 
(where x5, m (tan 2 ^) is minimized in Am 2 ) as a fuction 
of tan 2 9. From the figure it is possible to extract the 
maximum allowed values of tan 2 9 at any other CL. For 
instance, taking the global analysis in the LMA region 
from the figure we see that tan 2 6 = 3 (9 = 60°) is pos- 
sible only for A% 2 > 16.6 which corresponds to 99.98 for 
2 dof ( 3.65cr ). We present our results according to two 
different statistical criteria used in the literature in the 
defintion of the allowed paramers. The use of each crite- 
rion depends on the physics scenario to which the result 
of the analysis is to be applied. 

• Criterion 1 (CI): The regions at certain CL are 
defined in terms of shifts of the \ 2 function for 2- 
d.o.f, Ax 2 =4.6 (9.2) at 90 (99) % CL, with respect 
to the global minimum in the full plane. This cri- 
terion is used, for instance, in Refs. [1,6,18,9,8,19]. 
It is applicable when no region of the parameter 
space SMA, LMA, or LOW is a priory assumed to 
be the right one. 

• Criterion 2 (C2) : The regions at certain CL are de- 
fined in terms of shifts of the x 2 funcion for 2-d.o.f 
with respect to the JocaJ minimum in the corre- 
sponding region. This criterion holds when assum- 
ing that a given solution, SMA, LMA, or LOW, is 
the valid one. It clearly yields less restrictive limits. 
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FIG. 2. Xmin(tan 2 S) — xLn in the LMA (full lines) and 
LOW (dashed lines) regions as a fuction of tan 2 8, coming 
from the fit to the total rates (upper curves) and from the 
global data analysis (lower curves). 

This criterion is used, for instance, in Refs. [20,8]. 

In the figures we mark with a star the location of the 
minima used to define the contours in each case. For 
instance, in Fig. 1 when using criterion I, the contours 
arc defined with respect to minimum Xmin = 0.37/1 dof 
which is obtained in the small mixing angle solution. For 
criterion 2 we define each region with respect to its local 
minimum whose value and the corresponding goodness 
of the fit is given in Tables I and II respectively. Notice 
that the relatively bad goodness of the fit for the LMA 
and LOW regions from the analysis of the rates only is 
the reason why for example in Fig. 1 the LMA and LOW 
regions arc quite different using CI and C2. 

As seen in Fig. 1 and Tables I and II when using crite- 
rion 1, we do not find any solution in the second octant 
at 99 % CL from the analysis of the rates only. One must 
increase the CL to 99.6 (99.6) for the LMA (LOW) region 
to extend into values tan 2 9 > 1. This is in agreement, 
for instance, with the results of Ref. [16]. In Ref. [1] they 
find some small allowed region in the second octant for 
the LOW solution at the 99 % CL. We have traced the 
origin of this small discrepancy to their use of the ex- 
ponential approximation for the electron number density 
in the sun. In our calculation we use the solar neutrino 
fluxes from Ref. [21] and the new numerical parametriza- 
tion of the sun density as given by Bahcall [22]. We 
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TABLE I. Xmin an d its gof(goodness of the fit), maximum 
allowed values of tan 2 (8) at 90 (99) % CL and the confidence 
level at which 8 > ^ is allowed in the LMA region for the 
different set of observables and the two statistical criteria dis- 
cussed in the text. 
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TABLE II. Xmin an d its gof(goodness of the fit), maximum 
allowed values of tan 2 (8) at 90 (99) % CL and the confidence 
level at which 8 > ^ is allowed in the LOW region for the 
different set of observables and the two statistical criteria dis- 
cussed in the text. 



have explicitely verified that if using the exponential ap- 
proximation in our calculation the allowed LOW region 
extends into the second octant at 99 % CL. 

The situation is changed when the regions are defined 
according to criterion 2, as seen in Fig. 1 and Tables I 
and II. In particular both the LMA and LOW regions 
overlap at 99% CL because the value of x 2 m between the 
two regions is below the 99 % CL defined respect to the 
local LOW miminum (for this reason the LOW region in 
the figure is only shown at the 90 % CL). Notice also that 
in our representation we have chosen to cut the parameter 
space at Am 2 > 10~ 8 eV 2 . Recently, in Ref. [23] it has 
been pointed out that matter effects may be relevant for 
lower values of Am 2 . One must notice, however that for 
such lower mass values the simple analytic expresions in 
Eqs.(3) and (4) start loosing validity [15,16,24]. 

In Figs. 3 and 4 we show the allowed regions when ei- 
ther the data on the zenith angular dependence or the 
recoil electron energy spectrum are combined with the 
results from the total rates. The corresponding values 
of the absolute minimum of the x 2 fuction for the com- 
bination of rates plus zenith angular dependence data 
(rates plus recoil electron energy spectrum) are Xmin = 
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5.9/3 dof (22.1/15 dof) which are obtained for the SMA 
(LMA) solutions and are used in the construction of the 
allowed regions for criterion 1. In the figures we also 
show the corresponding excluded regions at 99 %CL by 
the new observables. 

As seen in Fig. 3 and in the tables the inclusion of 
the data on the zenith angular dependence of the Super- 
Kamiokande events leads, in general, to an increase in the 
maximum allowed values of tan 2 9 and a better CL for 
the second octant for both LMA and LOW regions. This 
is due to the fact that the best fit point for the angular 
distribution is obtained in the second octant (Am^ = 
3.7 x 1(T 6 eV 2 , tan 2 (6») = 5.9 with X m in,zen = L5 ) as 
shown in Fig. 3 [8]. 

This fact has been obviated in past analysis in the two 
neutrino oscillations and it must be taken into account 
to do properly the \ 2 analysis when the data on zenith 
dependence is included. In this respect, one must no- 
tice, for instance, that in their preliminary analysis on 
the zenith angle dependence [6] SuperKamiokande ob- 
tains the minimum \ 2 a t the boundary of their "cut" 
mixing parameter space sin 2 (29) = 1. In this way, their 
Xmin i s higher than the "true" minimum which is missed. 
Although, at present, the difference in the excluded re- 
gion defined with respect the "true" or with respect to 
the "cut" minimum is small, this can be of further im- 
portance when more data is acummulated. 

We have also explicitely verified that the results on the 
zenith angle excluded region as well as the best fit point 
position are very mildly dependent on the exact profile of 
the Earth. Very similar results can be obtained by using 
the analytical expressions valid for the two-step Earth 
density profile [10]. 

The effect of the inclusion of the recoil electron en- 
ergy spectrum is shown in Fig. 4. In this case when 
using criterion 1 we also find an increase in the maxi- 
mum allowed values of tan 2 9 and a better CL for the 
second octant for both LMA and LOW regions as com- 
pared to the results for the analysis of the rates only. 
However this is not the case when using criterion 2 as 
can be seen in the tables. The allowed regions from the 
global analysis are displayed in Fig. 5. We see from the 
figure that at 99 % CL both the LMA and the LOW re- 
gions extend into the second octant when using any of 
the two statistical criteria. In Fig. 2 we show the value 
of Xmin (tan 2 9) — Xmin m tne LMA and LOW regions 
(where x^^tan 2 9) is minimized in Am 2 ) as a fuction 
of tan 2 9 coming from the fit to the total rates (for which 
we recall that Xmin = 0.37) and to the global data set 

(for which Xmin = X 2 min ,LMA = 27 - )- As seen in the fig- 
ure the inclusion of the new observables leads to a small 
shift in the position of the local best fit points towards 
slightly larger values of tan 2 9 . However, as discussed 
before this not always translates into an increase on the 
maximum allowed value of tan 2 9 as the minimum also 
becomes deeper. 
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FIG. 3. Same as Fig. 1 but including also the data on 
the zenith angle distribution observed in Super-Kamiokande. 
The shadowed area represents the region excluded by the 
zenith angle distribution data at 99% CL. The location of 
minimum of the x 2 function for the zenith angle data is mar- 
ket with a point. 
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To summarize, in this paper we have studied the exten- 
sion of the second octant of the solar neutrino parameter 
space after including in the analysis all measured total 
event rates as well as all the 825-day Super-Kamiokande 
data on the zenith angle dependence and the recoil elec- 
tron energy spectrum of the events. We also have discuss 
the dependence of the results on the statistical criterion 
used in the construction of the allowed regions. Our re- 
sults are summarized in Tables I and II where we show 
the maximum allowed values of the tan at 90 and 99% 
CL when the different observables are included as well as 
the CL for the second octant of the parameter space. We 
have shown that when the zenith angle and the spectrum 
information is included the regions extend slightly more 
into the second octant. For the LMA this behaviour is 
mainly driven by the data on the zenith angular depen- 
dence since the best fit point for the zenith angle distri- 
bution is in the second octant. 

Finally just to comment that, the existence of solu- 
tions to the solar neutrino problem for 6 > j is not only 
of academic interest to the extent that in general they 
are perfectly allowed by the models of neutrino masses 
[25]. A particular interesting recent example of a predic- 
tive model can be found in Ref. [26] where, despite the 
predictivity of the model, solutions in both octants are 
equally possible. Furthermore, in the context of three 
neutrino mixing, the determination of the sign of Am 2 
of the LMA solution can be of relevance in the exact 
determination of the CP-violating phase [27]. 
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